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We study the dynamical properties of a two-dimensional ensemble of self-propelled dumbbells 
with only repulsive interactions. This model undergoes a phase transition between a homogeneous 
and a segregated phase and we focus on the former. We analyse the translational and rotational 
mean square displacements in terms of the Peclet number, describing the relative role of active 
forces and thermal fluctuations, and of particle density. We find that the four distinct regimes of 
the translational mean square displacement of the single active dumbbell survive at finite density 
for parameters that lead to a separation of time-scales. We establish the Peclet number and density 
dependence of the diffusion constant in the last diffusive regime. We prove that the ratio between 
the diffusion constant and its value for the single dumbbell depends on temperature and active force 
only through the Peclet number at all densities explored. We also study the rotational mean square 
displacement proving the existence of a rich behavior with intermediate regimes only appearing 
at finite density. The ratio of the rotational late-time diffusion constant and its vanishing density 
limit depends on the Peclet number and density only. At low Peclet number it is a monotonically 
decreasing function of density. At high Peclet number it first increases to reach a maximum and 
next decreases as a function of density. We interpret the latter result advocating the presence of 
large-scale fluctuations close to the transition, at large enough density, that favour coherent rotation 
inhibiting, however, rotational motion for even larger packing fractions. 

PACS numbers: 05.70.Ln, 47.63.Gd, 66.10.C- 


I. INTRODUCTION 

Active matter includes different kinds of self-driven 
systems which live, or function, far from thermodynamic 
equilibrium, by continuously converting internal energy 
sources into work or movement [1,2]. Nature offers many 
examples of this kind of condensed matter, at very dif¬ 
ferent scales: the cytoskeleton, bacterial colonies and al¬ 
gae suspensions, bird flocks and schools of fish are just 
some among others fH-flol|. Self-propelled units can also 
be artificially realized in the laboratory in many differ¬ 
ent ways, for example, by surface treatment of colloidal 
particles [Tlj . 

Active matter is inherently out of equilibrium and ex¬ 
hibits non-trivial properties that have no analogue in 
passive, equilibrium materials. For example, large scale 
coherent motion and self-organised dynamic structures 
have been observed in colonies of bacteria in the absence 
of any attractive interaction [T3 [pj . In addition, a phase 
separation into an aggregate and a gas- like phase has 
been found in theoretical models [19| - [27| and, recently, 
also in experiments [28| on suspensions of self-motile par¬ 
ticles only subject to steric interactions. 
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The dynamical properties of a suspension are signifi¬ 
cantly affected by self-propulsion. For instance, the com¬ 
parison of the diffusion constant D of a dilute solution of 
passive spherical colloids with the one of run-and-tumble 
bacteria shows that the Stokes-Einstein formula and the 
fluctuation-dissipation theorem do not hold for the active 
system mm- Indeed, assuming the Stokes-Einstein 
relation, D is given by D = ksT / (3irr)cr) where T is 
the temperature, rj the fluid viscosity, and a the diam¬ 
eter of the colloids. The diffusion coefficient in a three- 
dimensional system of run-and-tumble bacteria @ is eval¬ 
uated as Da ~ i„„/(6r) where r is the duration of each 
run and l run its length. Using t ft 1 s, l run ft 10 — 30 /im, 
water viscosity ry ft 10~ 3 Pa s, and cr ft 1 /im, one finds 
D a /D ft 10 2 . 

This simple order-of-magnitude argument suggests to 
analyse more carefully how the diffusive behaviour is af¬ 
fected by self-propulsion. Actually, several experimental 
studies addressed this question. Wu and Libchaber con¬ 
sidered the mean-square displacement of passive tracers 
coupled to a dilute suspension of Escherichia coli and 
found that an initial super-diffusive behaviour crosses 
over to normal diffusion at late times [Tsij]. The super- 
diffusive behaviour was interpreted as being due to the 
presence of coherent structures in the bacterial bath. A 
similar study was carried on by Leptos et al. on a sus¬ 
pension of algae. These authors found a linear time- 
dependence of the passive tracer’s mean-square displace¬ 
ment at all measured time-lags [30]. The tracer’s diffu- 
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sive constant was found to depend linearly on the density 
of swimmers (30) while, reducing the dimensionality to a 
film-like geometry, the density-dependence was enhanced 
to the power 3/2 [3l|. Results on the diffusion coefficient 
of tracers in contact with bacterial suspensions with hy¬ 
drodynamics playing a relevant role were given in [32 h34| . 

Other studies focused on the mean-square displace¬ 
ment of the active particles themselves. A linear de¬ 
pendence of the diffusion constant of the active swim¬ 
mers on the so-called active flux (active swimmers den¬ 
sity times their mean velocity) was found for different 
kinds of swimmers in contact with a solid wall [35|. The 
simulations of Hernandez-Ortiz et al. show ballistic be¬ 
haviour crossing over to normal diffusion for the swimmer 
and passive tracer particles though with different density- 
dependence of the diffusion constant in the latter regime. 
The swimmer’s diffusion constant decreases with den¬ 
sity at low swimmer density while it increases at large 
values; the diffusion constant of the passive tracer has, 
instead, a monotonic dependence with density and it con¬ 
sistently increases with it [15] . The Lattice Boltzmann 
study in (3f| also shows a cross-over between ballistic and 
diffusive behaviour at relative high self-propelled parti¬ 
cle density and a super-diffusive regime associated to the 
formation of large scale clusters at low densities. 

A detailed study of diffusion properties in models for 
self-propelled particles at different densities for the in¬ 
teresting cases where self-propulsion has been recognised 
to be an interaction capable of producing a phase transi¬ 
tion is lacking in the literature. Some results for the case 
of self-propelled polar disks have been presented in [13] . 
The paper by Gregoire et al. [37|, who considered a model 
with particle interaction that favors alignment a la Vic- 
sek , can also be mentioned in this respect. These authors 
interpreted the super-diffusive behaviour of [l3] in terms 
of the cross-over found in their model close to the critical 
point. 

Swimmers typically have elongated shape. We have 
therefore decided to study the diffusive properties of ro¬ 
tational and translational degrees of freedom in a system 
of active dumbbells. This model was introduced in [38|] to 
describe the experimental behaviour of a bacterial bath 
coupled to colloidal tracers. A first study of its phase 
diagram appeared in iffil [2bt where it was shown that for 
certain (relatively high) densities and active forces the 
system phase separates into dense and loose spatial re¬ 
gions. A later work [39|] focused on the dumbbell effective 
temperature defined in terms of a fluctuation-dissipation 
relation [4Cj. However, a detailed analysis of the transla¬ 
tional and rotational mean-square displacements in the 
full range of time-delays and varying the three more rel¬ 
evant parameters (temperature, activity and density) si¬ 
multaneously, was not performed yet. We present such a 
complete analysis in the homogeneous phase here. 

The paper is organised as follows. In Section |TT] the 
dumbbell model is reviewed. The diffusion behavior of a 
single dumbbell is described in Section Hill In Section ITVl 
the numerical results for the interacting active system are 


presented. The phase diagram is analysed in terms of the 
Peclet number. This is a preliminary step needed to fix 
the region of the parameter space to be considered for 
the measurements of the diffusive properties. Then the 
translational and rotational mean-square displacements 
in the homogeneous phase are studied in all their dy¬ 
namic regimes. Special emphasis is put on the analysis 
of the parameter dependence of the diffusion coefficient 
in the asymptotic limit. A discussion will complete the 
paper in Section fVl 


II. THE MODEL 


A dumbbell is a diatomic molecule formed by two 
spherical colloids with diameter <7d and mass m d linked 
together. The atomic positions are noted iq and r 2 in a 
Cartesian system of coordinates fixed to the laboratory. 
The colloids are subject to internal and external forces. 

Typically, one assumes that there is an elastic link be¬ 
tween the colloids modeled by the finite extensible non¬ 
linear elastic force 


F fene — 


kr 


1 ~{r 2 /rl) 


(1) 


with k > 0. The denominator ensures that the spheres 
cannot go beyond the distance ro with r the distance 
between their centres of mass. An additional repulsive 
force is added to ensure that the two colloids do not 
overlap. This is the Weeks-Chandler-Anderson (WCA) 
potential ITlI 


Lwca {? ) 


F L j(r) - H LJ (r c ) r<r c 
0 r > r c 


with 


14 j (r) = 4e 



( 3 ) 


where e is an energy scale and r c is the minimum of the 
Lennard-Jones potential, r c = 2 1 / 6 ( jd. 

The active forces are polar. They act along the main 
molecular axis n, are constant in modulus pointing in the 
same direction for the two spheres belonging to the same 
molecule j42j, and read 

Fact = Fact n . (4) 


We take the interaction between the spheres in differ¬ 
ent dumbbells to be purely repulsive and of the same 
WCA form as for the two colloids composing one dumb¬ 
bell. 
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The dynamic equations for one dumbbell are 




m d ic. l+1 (t) 


-7* i(t) + Ffene(ri,i + l) + T], 
2N 

^ wca l ij | -p 
■T a 


\ - dV.lL r 


7T, ''b 


■ acti 5 


“7U+1 (t) - Ffene(ri,i+l) + T] l+l 

' F ac tj , 


yy r i+1> 


j=0 

i#*+i 


<9u + iq r i+ ij 


(5) 
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with * = 1,3, ...27V—1, ry = 17-iq, 77 j = |r*j| and V^ c& = 
Kvcawith V^vca defined in Eq. @. Once the active 
force is attached to a molecule a sense of back and forth 
atoms is attributed to them; F act is directed from the 
ith colloid (tail) to the i + 1th colloid (head). The active 
forces are applied to all molecules in the sample during all 
their dynamic evolution. F act changes direction together 
with the molecule’s rotation. 

The coupling to the thermal bath is modelled as usual, 
with a friction and a noise term added to the equation 
of motion. 7 is the friction coefficient and we do not 
distinguish friction along the main molecular axis and 
transverse to it, as done in some publications |43|. The 
noise r) is a Gaussian random variable with 


(Via(t)) = 0 , (7) 

(via{t)ri jb (t')) = 2"/k B TSijS ab S(t - t') , (8) 

with k B the Boltzmann constant and T the temperature 
of the equilibrium environment in which the dumbbells 
move, a and b label the coordinates in d dimensional 
space. An effective rotational motion is generated by the 
random torque due to the white noise acting indepen¬ 
dently on the two beads. We consider initial condi¬ 
tions at time t = 0 such that the initial angle 9q = 9{ 0), 
randomly distributed between [—7r,7r], has zero mean 
Nic = 0. 

The surface fraction is 


This parameter is also a measure of the ratio between the 
work done by the active force in translating the center 
of mass of the molecule by a distance of 2<7<j, and the 
thermal energy scale. Another important parameter is 
the active Reynolds number 


Re ac t 


^d-Fact 

<Td7 


( 11 ) 


defined in analogy with the usual hydrodynamic 
Reynolds number Re = Ly/v, where v is the kine¬ 
matic viscosity of a given fluid, representing the ratio 
between inertial and viscous forces. Here we set L —> a d , 
v -A Fact /7 and v ->■ 7 cr 2 Jm A . 


III. A SINGLE DUMBBELL 

Before studying the interacting problem with numer¬ 
ical simulations in Sec. IIVI we derive analytically the 
translational and rotational mean-square displacements 
of the single dumbbell. 

The equation of motion for the position of the centre 
of mass, r cm = (iq + r 2 )/2, of a single dumbbell is 

2 m d r cm (t) = - 271 - 

cm (t) + 2F act (t)+m ( 12 ) 

with the new noise £(f) = r; 1 (t) + r/ 2 {t) with vanishing 
average, (£ a (i)} = 0, and correlation 

= 4y k B T S ab S(t - t') . (13) 

This is the Langevin equation of a point-like particle with 
mass 2md, under a force 2F act , and in contact with a bath 
with friction coefficient 2y at temperature T. 

The equation of motion for the relative position of the 
two monomers, r = iq — r 2 , is 

m d r(t) = — 7 r(f) + 2F int (f) + C(*) (14) 


</> 



with the new noise C(t) = r] 1 (t) — » 7 2 (t) having zero av¬ 
erage, (Ca(t)) = 0 , and correlation 


with Sd the area occupied by an individual dumbbell, 
S the total area of the box where the dumbbells move 
and N their total number. The spring is supposed to 
be massless and void of surface. Therefore, in d = 2, 
Sd = 7 tct^/ 2 . We impose periodic boundary conditions 
on the two directions. 

The Peclet number, Pe, is a dimensionless ratio be¬ 
tween the advective transport rate and the diffusive 
transport rate. For particle flow one defines it as Pe = 
Lv/D , with L a typical length, v a typical velocity, and 
D a typical diffusion constant. We choose L —> a d , 
v —> F act /7 and D —> D)// = k B T/(2j) of the passive 
dumbbell to be derived below; then, 


2 odF act 

k B T 


( 10 ) 


(')) = 47 k B T S ab S(t - t') . (15) 

Note that the noises £ and £ are independent, 
= 0, for all a, b at all times. F int includes the 
elastic and repulsive forces internal to the single dumb¬ 
bell. 

Equation (TbO) controls the molecule’s elongation and 
its rotational motion while Eq. m determines the trans¬ 
lational properties of the dumbbell. The internal force 
Fj nt affects the elongation of the molecule while the ther¬ 
mal noise adds fluctuations to it but, more importantly, it 
applies an effective torque and induces rotations. Equa¬ 
tions m and m are coupled by the fact that F act acts 
along the axis of the molecule, the orientation of which 
changes in time in the presence of thermal fluctuations. 
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A. Elongation and rotation 

Let us call U|| the instantaneous unit vector pointing 
from monomer 1 to monomer 2 along the axis of the 
molecule, 9 the angle between U|| and an axis fixed to 
the laboratory, and u_l a unit vector that is perpendic¬ 
ular to U|| at all times. Using r = ru|| + run, with r 
the modulus of r, u|| = 6u± and u_l = —0U|| (note 
that we use here the Stratonovich discretisation scheme 
of stochastic differential equations [Hj] and we are thus 
entitled to apply the usual rules of calculus) [45[, one 
has 

m d (r - rO 2 ) = - 7 r + 2 F int + C|| , (16) 

m d (2f6 + rO) = --yrO + </± , (17) 

where we decomposed the noise into the parallel and per¬ 
pendicular directions, C = C||H|| + ClU_l. The relations 
between the unit vectors in the fixed laboratory and the 
co-moving frame are given by 


where all factors are evaluated at the same time and D = 
2 7 ksT. The explicit calculation yields 
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d(r) 

d£ 


7' 


m 

d t 


= 2 (F int ) + U 7 - 1 (r- 1 ) , 


= 0 . 


(25) 

(26) 


The first equation is independent of 9 but it involves the 
average of different functions of r. The second equation 
implies (9) = 9 0 = 9(0) and, as we will take random 
initial conditions with average [#o]ic = 0 , then [( 0 )]i C = 
0. The statistics of 9 can be further analysed from the 
equation for the angular variance 


7 A(0 2 ) =2 ^ 7 - 1 (r- 2 ). (27) 


If we assume that r does not fluctuate around cr d , other¬ 
wise stated, that the molecule is approximately rigid (48j . 


u x = cos 9 U|| — sin# uj_ , 

u y = sin# U|| + cosd uj_ . (18) 

With this, for any noise we write 

C = (Cx cos 9 + (y sin 9) U|| 

T(-Cz sin0 + ( y cos0) uj_ , (19) 

and 

C|| = Cx cos 9 + Cy sin6> , 

C_l = -Cx sin 9 + ( y cos 9 . (20) 

The system of equations m-m for r and 9 cannot 
be solved exactly. We will assume that the internal and 
viscous forces are such that the inertial contributions (all 
terms proportional to m<j) can be neglected. We then 
have 


7 r = 2Tint + Cx cos 6 » + ( y sin 9 , (21) 

jr9 = — ( x sin0-|-^ y cos 9. ( 22 ) 

Putting together r and 9 into a vector y = (r, 9), this set 
of equations reads 

y a = h a [y] + g a p[y] (p ( 23 ) 

where the index /3 is a Cartesian one, £1 = and £2 = 
<^ y , and the index a yields yi = r and 2/2 = 9. The 
components of the vector h and the matrix g can be 
easily read from Eqs. m and (1^1) . In the last term the 
noise appears multiplying a function of the stochastic 
variable y. 

One can now average Eqs. ED and (l22l) over the 
Cartesian white noise by using the rules of Stratonovich 
stochastic calculus described in [3 [4(1 H?} : 

d 

(dapCp) = D(g v p g a p) ( 24 ) 


r « er d , 


this equation implies angular diffusion 
(0 2 ) = 9 2 + 2 D R t 

with the angular diffusion constant 

n _ D 2 k B T 
^ R o 9 2 

7 2 °d 70-a 


(28) 


(29) 


(30) 


The same technique can be used to compute all moments 
of the angular variable and thus show that it is Gaussian 
distributed within the rigid molecule approximation. 

It is interesting to compare our expression for the ro¬ 
tational diffusion constant, that for the sake of clarity 
we call T>^ umb in this paragraph, with that for self- 
propelled hard rods, as described in the Langevin ap¬ 
proach by Baskaran and Marchetti (43j| . In the model 
of Ref. [43 the rotational diffusion constant is given by 
Dft d = (kBTm )/(/ 7 ) where I is the moment of inertia 
calculated along the main axis of the rod and 7 is a fric¬ 
tion coefficient. If we identify the latter 7 with our fric¬ 
tion coefficient, and we take for the dumbbell the moment 
of inertia of two point-like particles with mass m = m d 
and diameter a d = 2 R, I = 2 mi? 2 , then our expression 
for the rotational diffusion coefficient coincides with the 
one in [43]. Consider now a rod of length l with the 
same aspect ratio (l = 4R, a = 2 R) and the same total 
mass 2m d of the dumbbell. One finds / = 2m d R 2 A with 
A « 2.94 H3, so that D dumb = AD™ d > D™ d . On the 
other hand, in the limit of a very long rod (l R) one 
has / = ml 2 /6 and £>^ umb = (D^ d l 2 )/(12R 2 ), so that 
£)duinb ^ as expected. Therefore the comparison 

with the model of Ref [43] suggests that the rotational 
diffusion coefficient of a dumbbell is always larger than 
the one of a rod in a suspension. 






5 


B. The center of mass 

We now focus on the statistical properties of the centre 
of mass position and velocity that depend upon the active 
force. One readily solves Eq. m 


In the limit t tj, 


(Ar c 2 m )W=4^t 




(38) 


/ v 0 m d \ v 0 m d — 2 _ t 
r C m (t) = r 0 +--e "m 

V 7/7 

+ ^~ fdt' [1 — ( * _t,) ][2F act (t') + £(t 4 )] ,(31) 

^7 Jo 

Vcm (t) = v 0 e _m d 4 

+ ;^~ f dt' e -^ (t - 4,) [2F act (t') + £(i 4 )] , (32) 

J o 

with r 0 = r cm ( 0 ) and v 0 = v cm ( 0 ). 

FromEq. (|32l) and thanks to [(cos0)] ic = [(sin0)] ic = 0 
one finds [(v cm )]j c = 0 and, after some long but straight¬ 
forward integrations, 

2m d [(t; 2 m )]j c = k B T + f a ^ , (33) 

7Hj + ta ) 


with 


. rn d 

ti = - , 

7 

~ 2k B T 


4D pd 


(34) 


in the long time limit, beyond tj. The time scales ti and 
t a are independent of the active force, and they are the 
usual inertial time and a characteristic time associated to 
rotational diffusion in the passive dumbbell, respectively. 
We also observe that in the passive limit Eq. (1351) reduces 
to the equipartition theorem for the kinetic energy of a 
point-like particle having the total mass of the dumbbell 
2md- As, typically, t a tj, one has 


2m d [(v 2 cmx )] ic ^k B T+^^. (35) 

With a similar calculation, starting now from Eq. m, 
we calculate the mean-square displacement (MSD) 

(Ar cm )(^) = [((r C m(i + io) — r cm(io)) )\ic (36) 

with [... ]j c the average over initial conditions at time 
t = 0. to is a sufficiently long time after preparation such 
that the stationary dynamics have been established and 
the mean-square displacement is therefore independent 
of to. Henceforth, t denotes time-delay. 

In the limit t <C tj, 


where 




k B T 

2 7 


(39) 


is the diffusion constant in the passive limit, F act = 0, 
see [39|| . This equation presents several time scales and 
limits. (Similar calculations for an active ellipsoid were 
presented in [49}.) For f/ < t < f a one finds 


<ArL>(i)=4i)£i + 



(40) 


that can still be split into the passive diffusive limit 
(Ar^ m ) = 4DP^ t for ti -C t < t*, and a ballistic regime 
(Ar^ m ) = (Fact/y ) 2 t 2 for t* < t -C t a , where the time 
scale t* is given by 

= 4F c P d 7 2 = 2k B T~i = (4_\ 2 _o|_ 

F 2 ct F 2 ct \PeJ 4 £>& 

= G7‘“' (41) 

Note that these two intermediate regimes do not exist if 
the parameters are such that t* < tj or t* > t a . They 
can also be easily confused with super-diffusion t a with 
1 < a < 2 if they are not well separated (ij ~ t* ~ t a ). 
See [39} and Fig. |4] below for more details. In the large 
Pe limit one has t* <C t a . In the last time-lag regime 
t t a , we recover normal diffusion, 


( Ar L)(i) = 4D a t , 

with the diffusion coefficient 
k B T 


D A {F act ,T,(j)= 0) = 


27 


1 1 / F act cr d \ ' 

2 ^ k B T ) 


(42) 


(43) 


In terms of the diffusion constant of the center of mass of 
the passive dumbbell and the Peclet number the above 
equation reads D A ( Pe, <j> = 0) = (1 + Pe 2 / 8 ). 

In the figures with numerical results for the finite den¬ 
sity problem shown in the next Section we include data 
for (f> = 0 that correspond to the single dumbbell limit. 


IV. FINITE DENSITY SYSTEMS 


(Ar 2 m )(f) = 2 ( , u 2 mx ) t 2 , (37) 

where ( v 2 mx ) = [(u 2 mx )] ic is the velocity given in 

Eq. (l33l) . The factor two is due to the sum over the 
two Cartesian directions. 


In this Section we present our numerical results. We fo¬ 
cus on three issues: the phase diagram, the translational 
diffusion properties and the rotational diffusion proper¬ 
ties. Details on the numerical method used for solving 
the dynamical equations m are given in [39} . We set 
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to d = (id = fc/j = e = 1 in proper physical units, and 
r 0 = 1.5, k = 30, 7 = 10, assuring over-damped motion 
and negligible dumbbell vibrations. Depending on the 
plots we used between 15000 and 20000 dumbbells in the 
simulations. 


A. The phase diagram 

Aspects of the phase diagram and the dynamics of this 
system were already established in J24|, |25[ . It was shown 
in these papers that at sufficiently low temperature and 
large active force the system phase separates into gas-like 
spatial regions and clusters of agglomerated dumbbells. 

The model has three important energy scales, e, ksT 
and T'act cr d- There is another energy scale related to the 
elastic constant k but, since we consider stiff dumbbells 
in this paper, we do not vary k. Keeping also e and the 
other parameters listed at the beginning of this section 
fixed, we will try to determine whether the phase diagram 
depends upon fc^T and F act <Td separately or only though 
their ratio, i.e. through the Peclet number Pe, as usually 
assumed [22| . The other free parameter to be varied is 
the global density of the sample, <f>. 

In the four panels in Fig. Q] we show the probability 
distribution function, p, of the local density, (f>i oc , for 
four values of the active force, F act = 0.05, 0.1, 0.5, 1. 
Each panel contains results for the same set of Peclet 
numbers Pe(T, F act ) obtained by different combinations 
of temperature and active force. The system has global 
packing fraction (j) = 0.5. We used the same operative 
definition of the local density as in jH], IHj] . We divided 
the full system in square plaquettes with linear size lOcrd 
that is much smaller than the linear size of the full sample 
and big enough to sample <f>i oc correctly. We improved 
the statistics by sampling over many different runs of the 
same kind of system. 

At low Pe the system is homogeneous and p is peaked 
around (f> i oc = (f> = 0.5. The critical Pe at which the sys¬ 
tem starts aggregating is approximately the same in all 
panels, Pe ~ 65. Around this value the density distribu¬ 
tion p not only becomes asymmetric but starts developing 
a second peak at (j)\ oc > 0.5 that characterises the dense 
phase in the system. Snapshots of typical configurations 
at Fact = 0.1 and four values of Pe are shown in Fig. 2. 
The location of the central peak at Pe less than the criti¬ 
cal value is independent of all parameters (apart from </>) 
while the location of the peak at cj)\ oc > 0.5 is situated 
at different values of </>i oc for different F act and the same 
Pe (compare the different panels in Fig. [I]). The reason 
for this is that the strength of the interactions between 
the dumbbells under different F act is different as F act /e 
varies with F act . A larger active force permits the dumb¬ 
bells to be more compact, while a lower one favors looser 
clusters. 

We repeated the analysis above for the cases with total 


packing fractions (f) = 0.3, 0.4, 0.6. We found the same 
values for the densities of the separated phases at (j) = 
0.4, 0.6 and Pe = 100, 200, and at <j> = 0.3 and Pe 
= 200. At <j> = 0.3 and Pe < 100 the effects of the 
presence of the spinodal line require a more elaborate 
analysis of the phase diagram, as discussed in [25]. In 
Table H we report these density values for the cases with 
Fact = 1, 0.5, 0.1, 0.05 and Pe = 100, 200. As observed, 
the coexistence values get closer for smaller active forces 
even though the Peclet number remains the same. 


Pe 

Fact = 0.05 

O 

II 

+3 

o 
r a 

fen 

Fact = 0.5 

-^act — 1 

200 

0.37 

0.34 

0.21 

0.049 


0.70 

0.71 

0.80 

0.890 

100 

0.44 

0.41 

0.37 

0.096 


0.66 

0.68 

0.77 

0.870 


TABLE I. Density values of the two coexisting phases mea¬ 
sured from the histograms in Fig. [J] at tj> = 0.5. The first and 
the third lines refer to the dilute phase while the two other 
lines correspond to the aggregated phase. Similar values are 
obtained at <j> = 0.4, 0.6 (Pe = 100, 200) and at <j> = 0.3 (Pe 
= 200 ). 

In Fig. [2] we show four snapshots of the system config¬ 
uration. The active force is F act = 0.1 in all panels and 
temperature is increased from left to right and from top 
to bottom. The configuration in the upper-left panel (Pe 
= 200 ) shows phase separation with large scale clusters 
while the configuration in the lower-right panel (Pe = 
20) is clearly homogeneous. The case Pe = 100 is in the 
segregated phase while the one for Pe = 66 is close to 
critical. 


B. Translational diffusion properties 

In Ref. [l3] the diffusion properties of a tracer im¬ 
mersed in a bacterial bath were monitored. A cross-over 
between a super-diffusive regime at short time-delays and 
a diffusive regime at long time-delays was reported. The 
cross-over time was found to increase linearly with the 
density of the active medium, showing that the cross¬ 
over is not due to the tracer’s inertia but to the dynami¬ 
cal properties of the bacterial bath. We explore here the 
same issues by focusing on the MSD of the center of mass 
of the dumbbells, defined in Eq. (13^1) . We will consider, 
for the rest of the paper, sufficiently low Peclet numbers 
such that the system will always be in the homogenous 
phase even though fluctuation effects can be relevant, as 
we will see. 
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FIG. 1. Probability distributions for the local density (j>\ oc at Peclet numbers given in the keys in the form Pe(T, F ac t) for 
different temperatures, and Fact = 0.05, 0.1, 0.5, 1 in the different panels. The global density of the system is cf> = 0.5. 


1. Dumbbell trajectories 


2. Four dynamic regimes 


Several single dumbbell trajectories are shown in Fig. [3] 
for different values of the temperature and global den¬ 
sity, under the same active force .F ac t =0.1. The trajec¬ 
tories correspond to a total time interval that includes 
the late diffusive regime (see below). At low tempera¬ 
ture and global density (T = 0.005 and cj> = 0.1, up¬ 
per left panel) we see periods of long directional motion. 
These are reduced at higher global density (</> = 0.4 up¬ 
per right panel). Increasing temperature at <j> = 0.1 
(T = 0.05 and T = 0.5 lower left and right panels, re¬ 
spectively) the trajectories become more similar to the 
typical ones of passive diffusion. While the trajectories 
are very stretched at T = 0.005, they become the most 
compact in the intermediate case at T = 0.05 and again 
quite stretched in the last case at T = 0.5. This be¬ 
haviour corresponds to the non monotonic behavior of 
the translational diffusion constant of Eq. fl3l) in terms 
of temperature. It decreases going from T = 0.005 to 
T = 0.05 while it increases going from T = 0.05 to 
T = 0.5. The single dumbbell diffusion coefficient, as 
calculated from Eq. (H3l) .is Da = 0.050, 0.0075, 0.025 
for the cases at T = 0.005, 0.05, 0.5, respectively. 


In Fig. [4] we show the center of mass MSD normalised 
by time-delay in such a way that normal diffusion ap¬ 
pears as a plateau. The four panels display data at four 
temperatures, T = 0.005, 0.01, 0.05, 0.1, all under the 
same active force F act = 0.1. Each panel has five curves 
in it, corresponding to five different densities given in the 
key. In all cases m d = 1 and 7 = 10 implying t/ = 0.1. 
The characteristic times fj, t*, t a are shown with small 
vertical arrows in each panel. These plots show several 
interesting features: 

- In all cases there is a first ballistic regime (the dashed 
segment close to the data is a guide-to-the-eye) with a 
pre-factor that is independent of <f> and increases with 
temperature as given by Eq. (1371) (The case t -C f / of the 
single dumbbell.) 

- Next, the dynamics slow down and, depending on T 
and (f), the normalised mean-square displacement may 
attain a plateau associated to normal diffusion (low T ) 
or even decrease, suggesting sub-diffusion. (The case 
i/ <C f -C t* -C f a of the single dumbbell.) 

- The dynamics accelerate next, with a second super- 
diffusive regime in which the curves for all <f> in each panel 
look approximately parallel and very close to ballistic at 
T = 0.005, 0.01, 0.05. (The case f/ -C t* -C t -C t a of 
the single dumbbell.) 

- Finally, the late normal diffusive regime is reached with 
all curves saturating at a plateau that yields the different 












FIG. 2. Snapshots of the system for F ac t = 0.1 (cfr. the upper right panel in Fig. [l]), and T = 0.001, 0.002, 0.003, 0.01 
corresponding to Pe = 200, 100, both phase separated, Pe = 66, close to critical, and Pe = 20, homogeneous (from left to right 
and from top to bottom). The global density is tj> = 0.5. 


Da coefficients. (The case t/ <C t*, t a <C t of the single 
dumbbell.) 

It is hard to assert whether the intermediate regime 
is super-diffusive or simply ballistic as the time-scales t* 
and t a are not sufficiently well separated (and not even 
ordered as t* < t a in the last panel). Moreover, in the 
last two panels (high T or low Pe) the diffusion-ballistic- 
diffusion regimes are mixed, due to the fact that the con¬ 
dition t* -C t a is no longer satisfied. The effective slope 
in the intermediate super-diffusive regime decreases when 
the density increases. 

A rather good fit of the finite density data in the limit 
Pe 1 and for time-delays such that t > t* is achieved 
by using the single dumbbell expression in Eq. (1551) 

<ArL>(t) = - 1 ~ e , DRt ) , (44) 

without the first term (negligible if Pe 1) and upgrad¬ 
ing the remaining parameters, D A and D to be density- 
dependent fitting parameters, as done in [HI Ho|. This 
is shown in Fig. [5] (left panel). For not that large values 
of Pe one could recover the remaining parameter and use 

instead (Ar 2 cm )(t) = 4D^t + 4D^(t - *"* ”"* ) with 

U R 

an additional fitting parameter. Figure [5] (right panel) 


also shows a good agreement between the values of D R 
found in these fits and the values of the rotation diffu¬ 
sion coefficient Dn(F ac t, T, cf>) coming from the late time- 
delay diffusive regime in the rotational MSD discussed in 

sec.' irm 

The cross-over time-delay between the last ballistic or 
super-diffusive, and the diffusive regimes seems quite <jh 
independent in the first two panels T = 0.005, 0.01 and 
it increases, though rather weakly, with </>, in the last 
two panels, T = 0.05, 0.1, see the inclined dashed line in 
the last panel that is also a guide-to-the-eye. This cross¬ 
over time-delay is the one that we could associate to the 
cross-over time between a superdiffusive regime and the 
last diffusive regime found in the experiment in [l3j] . The 
strongest effect of density is though on the first diffusive 
or sub-diffusive regime. 

In summary, no large qualitative change in the cen¬ 
ter of mass MSD behavior is observed in the range 
4> € [0,0.6]. There is just a natural slowing down of the 
dynamics with larger packing fractions that translates 
into a change from diffusive to sub-diffusive behavior in 
the second regime and a general decrease of the diffu¬ 
sion constant in the last regime for all Pe. We study the 
dependence of the diffusion constant with Pe in detail 
below. 
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T=0.005, (>=0.1 T=0.005, (>=0.4 



FIG. 3. Several trajectories of one bead in different dumbbells. The active force is Fact =0.1 and the evolution time interval is 
t = 10000 in all panels. The other parameters are T = 0.005 and <f> = 0.1, 0.4 (upper panels, left and right), and T = 0.05, 0.5 
and 0 = 0.1 (lower panels, left and right). The values of the single dumbbell diffusion coefficient, as calculated from Eq. (1431) . 
are Da = 0.050, 0.0075, 0.025 for the cases at T = 0.005, 0.05, 0.5, respectively. 


3. The late-epochs translation diffusion coefficient 

Let us now discuss the normal diffusive regime at 
longest time-lags. In [39] we studied the translational 
diffusion coefficient Da as a function of F ac t and (j) at 
fixed temperature. In particular, we compared the <f> 
dependence to the Tokuyama-Oppenheim law for col¬ 
loids [50]. Here, we first examine, instead, the T and 
<fi dependence of Da for fixed active force, F act . Then 
we consider how the dependence of Da from T and F act 
can be re-expressed in terms of the Peclet number. The 
main results for Da obtained in [39| will also be revisited 
in this subsection. 

The first question we want to answer is whether Da de¬ 
pends on kgT as for the single dumbbell case ( (j) = 0), the 
functional form recalled in Eq. (|43[) . For T’actO’d/fesT 
1 such that the quadratic term can be neglected this 
equation implies the linear growth of Da with ksT as 
in the passive limit. Instead, when the second term dom¬ 
inates, i.e. for very small thermal energy with respect to 
the work performed by the active force, Da should decay 
as 1 /(fcgT) with a slope that is quadratic in F act <Jd- 

In Fig. [S] we display Da as a function of T for various 
values of (f> given in the key and F act = 0.1. The theo¬ 
retical values for cf> = 0 are included in the figure (with 
open triangles joined by a dotted curve). Here, we used 
the measured value for the distance between the centres 


of the two colloids, that is r « 0.96 era- In the rest of 
this section we simply call <Td the molecular length and 
we take era = 0.96. 


The error-bars are smaller than the symbol size and 
we do not display them. The curves show a minimum 
located at fcsT m ; n = F act ad/V2 for 0 = 0, that weakly 
increases with (f>. The two regimes, Pe <Cl and Pe 1, 
still exist and Da is dominated by thermal fluctuations in 
the former and by the work done by the active force in the 
latter as in the single dumbbell limit. We see a saturation 
of Da at small values of T for 0 > 0.2 and therefore 
the breakdown of the single dumbbell 1 /(ksT) behaviour 
at low temperatures. Instead, at high temperatures Da 
seems to retain the linear growth with temperature of the 
single dumbbell at least for the temperatures used in the 
simulations. 

Figure [HI also shows that for the Pe that we used Da 
is a decreasing function of (f> at all fixed temperatures. 
This fact can be better appreciated in the left panel in 
Fig. 0 where Da is plotted as a function of (j> for various 
temperatures given in the key. (Recall that the (f> depen¬ 
dence of Da at fixed T and for different active forces was 
discussed in 139] where it was shown how the Tokuyama- 
Oppenheim [501 ] law of the passive system was simplified 
under activation to a decay that is close to a simple ex- 
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FIG. 4. The center of mass MSD normalised by time-delay, for an active system under Fact = 0.1 at T = 0.005, 0.01, 0.05, 0.1 
(Pe = 40, 20, 4, 2), at the different densities given in the key. The Peclet number decreases from one panel to the other. It 
induces a strong qualitative change in (Ar^), see the text for a detailed discussion. The two dashes in the first panel represent 
the ballistic dependence ~ t 2 . The dashed line in the last panel is a guide-to-the-eye for the density dependence of the last 
cross-over time-delay that increases weakly with (j>. The horizontal dotted lines at long times correspond to the values of the 
single dumbbell diffusion constant D A from Eq. (S3). The vertical black arrows indicate the single dumbbell time-scales tj and 
t a , while the red arrows indicate the single dumbbell characteristic time t *, for each case. In the case with Pe = 4 the times t* 
and t a coincide. The curves for (j> = 0 are obtained from Eq. (1381) . valid in the limit t ti = 0.1; for this reason they start 
from the middle of the graph. 


ponential. We will come back to this issue below.) 

The non-monotonicity of Da as a function of T already 
discussed in Fig. [S] is confirmed by the data presentation 
in Fig. [3 with the minimum situated around T ~ 0.07. 
In the right panel we observe the opposite behaviour in 
the ratio D A (F act = 0.1,T,(j))/D A (F ac t = 0.1,T,0), first 
growing for increasing T to reverse its trend at around 
T ~ 0.05 — 0.07. Consistently with the behaviour found 
in HH, there are temperatures such that the data for 
the above ratio cross each other when the density is in¬ 
creased, see for example T = 0.01, 0.1 (or Pe = 2, 20). 
The right panel in Fig. [3 also shows that a very small 
density can have relevant effects on the behaviour of the 
diffusion coefficient. 

We have repeated this analysis for a stronger active 
force and we found that the results are consistent, with 
a cross-over temperature that grows with Factor, as pre¬ 
dicted by the single dumbbell equation, though we cannot 
assert that the dependence be linear. 

Next, we analyse in Fig. [8] whether the ratio of dif¬ 
fusion coefficients of the active system at finite density 


and single passive dumbbell F A (F act , T, <p)/Df.^ depends 
only on the Peclet number, as it does for the single dumb¬ 
bell problem. With this aim, we fix F act and we vary T, 
and the values Pe = 4, 20, 40, 66 in each panel are ob¬ 
tained from three different combinations of F act and T. 
In all panels the collapse of data is very good. Note the 
change in concavity of the collapsed data that occurs at 
Pe = 20. This value is relatively far from the transition 
between homogeneous and segregated phases estimated 
in [25], [3f|, and the system configurations are still ho¬ 
mogeneous, see the last panel in Fig. [3 though with a 
distribution of local densities, <f> i oc , with a certain width, 
see Fig. [3 

These results suggest 

D A {F acU T, 4 ) = k B T f A (Pe, 4) (45) 

with / A (Pe, 0) = ( 27) _1 (1 + Pe 2 /8) = 

F A (F act ,T, 0)/(fcsT) and f A a decreasing non-linear 
function of <f> at fixed Pe. This relation is equivalent to 

D A (F act ,T,4) = / A (Pe, 4) (/ , 6) 

D A (F act ,T, 0) / A (Pe,0) ' 1 ’ 
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FIG. 5. Left panel: Fit of the center of mass MSD normalised by time-delay with the expression in Eq. (l44l) . Pe = 40 and 
the different data sets correspond to the global densities given in the key. Right panel: The rotation diffusion coefficient as 
extracted from the fit of the center-of-mass MSD shown in the left panel (D^ in Eq. (1441) - blue stars), and from the late-time 
diffusive rotational MSD directly measured as shown in Sec. IIV Cl ( green squares). Statistical errors are within the size of 
the symbols used in the figure. 



FIG. 6. The center of mass diffusion coefficient, Da, as a function of temperature, T, at different densities <j> given in the key 
and fixed active force, Fact = 0.1. The triangular data points joined by a dotted curve labeled 4> = 0 are given by the theoretical 
Da with era = 0.96. 


The l.h.s. is what we studied in |39j as a function of 
.Fact and <j>, keeping T fixed, and we proposed 

-P. 4 (^acti T,(f>) _ -h(f act )0 (47'1 

D A (F act ,T, 0) 1 j 

with b a non-monotonic fitting function of F act . Knowing 
now that D A /(ksT) depends on F ac t and T only through 
Pe, we deduce 

D A {F act ,T, 0) = DU (1 + Pe 2 / 8 ) e~ b ^ . (48) 

Note that in [39j the maximum in b appeared at F ac t ~ 
0.1 that, for the temperature used, T = 0.05, corresponds 
to Pe ~ 4. Thus, f A ( Pe, 0) should be monotonically in¬ 
creasing with Pe, at all fixed 0, as it results when compar¬ 
ing the data on the different panels in Fig. [ 8 ] In Fig. [ 8 ] we 
included, with dotted black lines, the exponential fits in 


Eq. (H51) where the only free parameter is 5(Pe). The val¬ 
ues of 6 (Pe) are 1.1, 1.6, 2.8, 4.1 for Pe = 4, 20, 40, 66 , 
in agreement with what we reported in ]39|. 

However, while we see that the exponential fit is very 
good at all 0 for Pe = 40 and Pe = 66 , it is not as good 
for the smaller Pe data. The red line-points in Fig. [ 8 ] 
represent, instead, the result of the fit 

D A (F act , T, (j>) = D A (Fact , T, 0) 

x [1 + ai(Pe) (j) + a 2 (Pe) 0 2 ] . (49) 

This functional form gives a better representation of the 
data than the exponential for Pe = 4 and Pe = 20, which 
is, in a sense, natural since one expects to recover a rather 
complex Tokuyama-Oppenheim like form in the limit Pe 
—> 0. The exponential and polynomial fits are of equiva¬ 
lent quality for Pe = 40, while the polynomial fit is clearly 
worse than the exponential one for Pe = 66 . The fitting 
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FIG. 7. Center of mass diffusion constant at various T given in 
the one for one dumbbell (right panel) both as a function of <j>. 
and (f> = 0.001 in the right panel. The active force is fixed to F t 



the key (left panel) and center of mass diffusion constant over 
The smallest densities considered are (j> = 0 in the left panel 
, c t =0.1. The symbol convention is the same in the two panels. 






FIG. 8. Center of mass diffusion constant, Da, over fesT for Pe = 4, 20, 40, 66 as shown in the keys. For each value of the 
Peclet number three couples of values of F ac t and T were used. The black dotted lines are exponential fits to the data points, 
as suggested by Eq. (|48|. The red line-points are the quadratic fit in Eq. (1491) with fitting parameters given in the keys. See 
the text for a discussion. 


parameters are given in the keys. One notices that ai is C. Rotational diffusion properties 

negative in all cases while a 2 changes sign from negative 

at Pe < 20 to positive at Pe > 20 (leading to a growing 

behaviour at large <f> that is not physical). At Pe = 20 

the density dependence is almost linear as 02 is very close 

to zero. 


Having discussed in detail the translational diffusion 
properties we turn now to the rotational ones. 
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FIG. 9. The angular MSD for the Pe numbers in the labels and, with different lines, various densities given in the keys to each 
panel. The dash in the first panel highlights the initial ballistic behavior. The vertical arrows indicate the characteristic times 
ti, t *, t a ■ The horizontal dotted lines at long times correspond to the values of the single dumbbell diffusion constant Dr 
from Eq. flUD - Note that the vertical scale is different in the two panels above and in the two panels below. 


1. Dynamic regimes 


In Fig. [9] we display the angular MSD normalized 
by time-delay. The four panels show data obtained for 
the same parameters as the ones used in Fig. 0] with 
F a ct = 0.1. Each panel, corresponding to the cases 
with T = 0.005, 0.01, 0.05, 1 (Pe = 40, 20, 4, 2, 
respectively), includes curves for five finite densities, 
(j> = 0.1, 0.2, 0.3, 0.4, 0.5, and the single dumbbell 
limit, cj) = 0, as labeled in the key. These plots also show 
several interesting features: 

In all cases there is a first ballistic regime with a pre¬ 
factor that is independent of <j> and increases with temper¬ 
ature (The case t -C tj = m &/7 of the single dumbbell.) 
- Next, the dynamics slow down and, depending on T 
and the normalised MSD may attain an ever-lasting 
plateau associated to normal diffusion for low cj) at any 
temperature, or even decrease, suggesting sub-diffusion, 
at high enough </>. 

At low temperature T = 0.005, 0.01 (Pe = 40, 20) and 
sufficiently high density the dynamics accelerate next, 
with a second super-diffusive regime that crosses over to 
a final diffusive regime. 

In the late normal diffusive regime all curves saturate 


and the height of the plateau yields the different Dr co¬ 
efficients that we discuss below. 

The effect of Pe and (f> are stronger on the rotational 
MSD than on the translational MSD. New regimes ap¬ 
pear in the rotational collective motion with respect to 
the individual molecular limit. In the phase separated 
regime the dumbbell clusters rotate [24], [25( • It is possi¬ 
ble that strong fluctuations not far from the critical point 
(Pe = 20, 40) have an important rotational component 
than enhances/advects rotational diffusion giving rise to 
an observable contribution to displacement also manifes- 
tating itself in the appearing of new dynamical regimes. 


2. The late-epochs rotation diffusion coefficient 

We now study whether the linear temperature depen¬ 
dence of the single dumbbell angular diffusion constant, 
Eq. (1301) . survives the interactions between dumbbells in 
the finite density problem, see Fig. [10] The data points 
are compatible with a linear behaviour at sufficiently 
high temperature, with a slope that depends upon <f>. 
The trend in the curves reverses below the cross-over at 
T ~ 0.01 with larger values of Dr for larger values of </> 
(see the right panel in the same figure). 
















14 




FIG. 10. Angular diffusion constant as a function of T, for F ac t = 0.1 at fixed 0 in logarithmic scale (left panel) and a zoom 
on the low temperature behaviour in linear scale (right panel). The data points joined with a line for 0 = 0 represent the 
theoretical expectation (30]). 


From Fig. [TT] one easily concludes that the Fact- 
independence of Dr is lost as soon as the interaction 
between dumbbells is switched on at finite density. This 
fact can be seen, for instance, by comparing the T = 0.1 
data, one of the two temperatures included in both pan¬ 
els, sharing the same value, slightly larger than 10 -2 , 
at 0 = 0. While in the case F act = 0.1 (left panel) 
Dr clearly decreases with 0, in the case F act = 1 (right 
panel) Dr is almost constant. These figures also show 
the change in trend operated at an F act -dependent T: 
at high temperature Dr decreases with 0 while at low 
temperature Dr increases with 0. The change occurs at 
T ~ 0.01 for F act = 0.1 and at T ~ 0.1 for F ac t = 1 
suggesting that the change is controlled by Pe. 

Finally, we analyse whether Dr/ (k B T) depends on 
F act only via Pe. To this end, in Fig. [12] we repeat the 
analysis shown in Fig. [8] for Da- The four panels show 
D R /(k B T) against 0 for Pe = 4, 20, 40, 66. In each 
panel we include data for three pairs of F ac t and T lead¬ 
ing to the same Pe. We see that the data points collapse 
on different master curves in each panel. This suggests 

F fl (F act ,T, 0) = k B T f R { Pe, 0) (50) 

with Jr( Pe, 0) = /r(0, 0) = 2 /(7 <t 2 ). The data also show 
a change in trend of the function / r at around Pe = 20. 
At low densities, while the master curve decreases with 
0 for Pe < 20, it becomes flat at Pe = 20 and it increases 
with 0 for Pe > 20. This would suggest: 

jA(Pe, 0) — + a(Pe, 0) , (51) 

700 

with a(Pe, 0) almost linear in 0 and the slope changing 
sign at Pe ~ 20 for small 0. All panels, i.e. at all Pe, 
show a cross-over at high enough densities after which 
the rotational diffusion constant decreases with increas¬ 
ing density. 

A possible explanation of the different density- 
dependence of Dr at small and large Peclet can be found 


from following the evolution of a single tracer dumbbell 
at intermediate densities, 0 r; 0.4 for example, as it can 
be seen in the supplementary movies in Ref. ;51]. One 
observes that at low Pe (Pe = 2) the system is very uni¬ 
formly distributed and the movement of the tracer dumb¬ 
bell is inhibited by the ‘cages’ formed by surrounding 
dumbbells. Collisions are frequent but each of them only 
produces a small angular displacement. In this case the 
effect of increasing the density is to decrease both the ro¬ 
tational and translational diffusion coefficients. On the 
other hand, at high Peclet (Pe = 40), small fluctuating 
clusters can be observed (their presence is also signalled 
by a peak in the structure factor )39|). This has rel¬ 
evant effects on the behaviour of the tracer dumbbell. 
First, there are particle depleted regions which are large 
enough to allow significant angular displacements with¬ 
out collisions. Second, angular displacements appear to 
be enhanced when the tracer dumbbell meets a cluster 
and is advected by its motion. On the other hand, at still 
higher densities the cage effect becomes again preeminent 
so that rotations are inhibited and Dr decreases. Note 
that both Dr and Da change behavior at Pe = 20 (the 
translational diffusion coefficient Da is a convex func¬ 
tion of density for Pe < 20 and changes curvature for Pe 
> 20). We find the fact that these cross-overs occur at 
the same Peclet worth to be stressed even though it is 
difficult to argue about its implications. 


V. CONCLUSIONS 

We presented a thorough study of the translational and 
rotational MSD of a system of interacting active dumb¬ 
bells. We focused on the regimes where the global system 
is homogeneous. Higher densities than the ones used 
in [T3 . [ 30 ] have been considered with the Peclet number 
small enough (possibly much smaller than in the experi¬ 
ments) to keep the system in the homogenous phase. 
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FIG. 11. Angular diffusion constant for various T under the same active force _F ac t =0.1 (left panel) and F ac t = 1 (right panel). 
The only common temperatures on the two panels are T = 0.05, 0.1 and, consistently, Dr at these temperatures is the same 
at <j) = 0 . 
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FIG. 12. Rotational diffusion constant over fcsT for Pe = 4, 20, 40, 66 as shown in the keys. For each value of the Peclet 
number three couples of values of F ac t and T were used. All master curves take the value 2 /( 7 er|) « 0.215 at <j> — 0. At small 
<j> , there is a noticeable change in trend at Pe ~ 20. At sufficiently large 1 p all curves decrease with increasing (j>. 


We first analysed the single molecule dynamics as a 
benchmark to later characterize the finite density ef¬ 
fects. In the passive case, Pe = 0, the translational 
and rotational MSDs show a standard cross-over from 
ballistic motion to normal diffusion at the inertial time 
tj = rrid/'Y- Under the active force, the normal diffu¬ 
sion of the center of mass is accelerated after a time-scale 
t* oc t a / Pe 2 with t a = 7(Jj/(2A:bT) and, still later, after 
t a , a new diffusive regime is reached with a diffusion con¬ 
stant that is enhanced with respect to the one in the pas¬ 


sive limit as a quadratic function of the Peclet number. 
Instead, the rotational properties of the active dumb¬ 
bell are not modified by the longitudinal active force; all 
torque is exerted by the thermal noise. 

Then we turned to the analysis of the mixed density 
and active force effects on the collective motion of the 
interacting system. 

The rich dynamic structure of the center of mass trans¬ 
lational motion of the single molecule, with the four dis¬ 
tinct time regimes summarized above, survives under 
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finite densities with modified parameters. The super- 
diffusive behaviour shown in [131 ] is reminiscent of the 
second ballistic regime in the interacting active dumbbell 
system at finite densities. The diffusion constant Da in 
the last diffusive regime has a non-monotonic dependence 
on temperature, as for the single dumbbell case, and it 
decreases with increasing self-propelled particle density 
at all temperatures. Moreover, the ratio Da/^bT) de¬ 
pends on temperature and active force only through the 
Peclet number at all densities explored. This ratio, at 
fixed density, is an increasing function of Pe. All these 
results are consistent with those found in our previous pa¬ 
per )39] where it was also shown that the ratio between 
the translational coefficient diffusion at finite density and 
the one for the single dumbbell had a non-monotonic Pe 
dependence. 

Next we moved to the analysis of the rotational MSD. 
While in the single dumbbell case its time-delay depen¬ 
dence is rather simple, with a single cross-over between 
ballistic and diffusive behaviour, intermediate regimes, 
roughly for ij <C i <C i* and t* -C t -C t a , appear 
at finite densities. The late epochs diffusion constant 
Dr increases with temperature (though not linearly) at 
all densities and active forces simulated. The indepen¬ 
dence on active force is lost at finite densities. The ra¬ 
tio Dr/{ ksT) depends on temperature and activity only 
through the Peclet number. At low densities, its depen¬ 
dence on density changes from decreasing at low Pe to 


increasing at high Pe. This change in behaviour can 
be related to the large scale density fluctuations that ap¬ 
pear close to the transition from the homogeneous to the 
aggregated phase at a critical Pe. In the aggregated 
phase lar ge a nd rather compact clusters rotate coher¬ 
ently [2J, |25[. Not far from the transition, in the ho- 
mogenoues phase, fluctuating clusters with some coher¬ 
ent rotation are observable and these may be the cause 
for the increase of Dr with </>. On the other hand, at large 
enough densities rotations are strongly inhibited and the 
value of Dr decreases for all Pe. 

The fluctuations of translational and rotational dis¬ 
placements have been characterized in [52] . Special em¬ 
phasis was put on the identification of the regimes in 
which the fluctuations are non-Gaussian. See this refer¬ 
ence for more details. 

After this work we plan to analyse the motion of trac¬ 
ers in contact with this active sample and, especially, to 
analyse the existence of a parameter to be interpreted as 
an effective temperature from the mobility and diffusive 
properties of the sample and the tracers, in the manner 
done in [29|, [53l46T| for different active systems. 
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